Hole-doped semiconductor nanowire on top of an s-wave superconductor: A new and 
experimentally accessible system for Majorana fermions 
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Majorana fermions were envisioned by E. Majorana in 1935 to describe neutrinos. Recently 
it has been shown that they can be realized even in a class of electron-doped semiconductors, on 
which ordinary s-wave superconductivity is proximity induced, provided the time reversal symmetry 
is broken by an external Zeeman field above a threshold. Here we show that in a hole-doped 
semiconductor nanowire the threshold Zeeman field for Majorana fermions can be very small for some 
magic values of the hole density. In contrast to the electron-doped systems, smaller Zeeman fields 
and much stronger spin-orbit coupling and effective mass of holes allow the hole-doped systems to 
support Majorana fermions in a parameter regime which is routinely realized in current experiments. 

PACS numbers: 74.78.-w, 03.67.Lx, 71.10.Pm, 74.45.+C 



Recently some exotic condensed matter systems, 
such as the Pfaffian states in fractional quantum 
Hall (FQH) systems chiral p-wave superconduc- 

tors/superfluids jH-lol, topological insulator (TI) [Toj , 
as well as a ferromagnet-superconductor heterostructure 
[lH, have been proposed as systems supporting quasi- 
particles with non-Abelian statistics [l2l]. These sys- 
tems allow a special type of quasiparticles called Ma- 
jorana fermions which involve no energy cost. The sec- 
ond quantized operators ji for the Majorana excitations 
are self-hermitian, = (particles are their own anti- 
particles), which lies at the heart of their non-Abelian 
statistical properties. Due to the fundamental difference 
of Majorana fermions from any other known quantum 
particles in nature, the emergence of these particles in 
solid state systems would in itself be an extraordinary 
phenomenon. Their potential use in fault-tolerant topo- 
logical quantum computation (TQC) [12] makes their re- 
alization in controllable solid state systems even more 
significant. 

It has been shown recently 13|-18| that an electron- 
doped semiconducting thin film or nanowire with a siz- 
able spin-orbit coupling can host, under suitable condi- 
tions, Majorana fermion excitations localized near de- 
fects. This proposal followed on an earlier similar pro- 
posal in the context of cold atomic systems pjj. When 
the film or the nanowire is in the presence of a Zeeman 
splitting Vz (with Lande factor ^*) and an s-wave su- 
perconducting pair potential A, which can be proximity 
induced by a nearby superconductor, the system enters 
into a topological superconducting (TS) state for 
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Here ji is the chemical potential in the semiconductor 
which is controlled by the density of doped electrons. 
Despite the theoretical success, the requirement Eq. ([T]) 
for the TS state in an electron-doped nanowire leads to 
two obvious experimental challenges: a low electron den- 
sity and a high magnetic field. For a small carrier density 



a nanowire tends to become insulating due to the strong 
disorder-induced fluctuations of the chemical potential. 
A high magnetic field, on the other hand, can be detri- 
mental to pairing as well as 5- wave proximity effect itself. 

In this Letter we show that a hole- doped semiconductor 
nanowire can solve all these problems encountered in the 
electron-doped systems. The hole-doped nanowire is very 
different in many respects from its electron-doped coun- 
terpart due to its different band structure and the value of 
the effective spin of the carriers. For some "magic" values 
of the carrier (hole) density the threshold Zeeman split- 
ting for the TS states and Majorana fermions can become 
very small, therefore the constraint on the carrier density 
as given in Eq. ([T]) is absent for the hole-doped nanowires. 
Furthermore, the effective mass and spin-orbit coupling 
in the p-type valence band holes are much larger than 
electrons, which leads to a larger Fermi vector kp. This 
larger kp leads to a larger required carrier density (~ 10^ 
cm~^) for the TS state, which, remarkably, is now rou- 
tinely achieved in many experiments |2QM22I]. The large 
carrier density provides strong screening of the disorder 
potentials, leading to much smaller fluctuations of the 
chemical potential [23| in the nanowire. Furthermore, the 
small ratio between the Zeeman coupling and the spin- 
orbit energy (orders of magnitude smaller than that in 
the electron-doped systems) leads to a small carrier mo- 
bility requirement for the hole-doped TS state (3 order 
of magnitude smaller than that for the electron-doped 
TS state), as pointed out recently in [2^. Let us also 
point out that the superconducting proximity effect on a 
hole-doped nanowire has been observed in recent experi- 
ments [22]. It seems therefore that a Majorana-carrying 
TS state is tantalizingly close to experimental reach in a 
hole-doped nanowire. 

Set-up and Hamiltonian: The experimental setup is 
illustrated in Fig. [TJi, where a hole-doped semiconductor 
nanowire is placed on top of an s-wave superconductor. 
We focus on large-band gap semiconductors, where the 
single particle Hamiltonian of holes is described by the 
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FIG. 1: (a) Schematic plot of the experimental setup, (b) 
Illustration of the valence band structure of ID hole-doped 
nanowire with a Zeeman field. Only the lowest band along 
the y and z directions is considered. There are two heavy 
hole (dashed black and solid green curves) and two light hole 
(blue dash-dotted and pink dotted curves) bands. The thick 
red line (amplified as shadow in inset) gives the regime of 
chemical potential with a single Fermi surface, which leads 
to the topological superconducting state. The parameters are 
for a hole-doped InAs nanowire with 71 = 20, 72 = 8.5. a = 
3.3 X 10^ m/s, Vz 1.5 meV, 14 nm, Ly 10 nm. 



four-band Luttinger model [25| (henceforth we set h ■ 
and m = —1), 
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7i , 



)V^--f2{V-Jf-ia{3xV)-z + V,J,-fi, 
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where a is the Rashba spin-orbit coupling due to the 
inversion- symmetry breaking, and the fourth term is the 
Zeeman field Vz = q^iibB generated by the external 
magnetic field along the z direction. J is the total an- 
gular momentum operator for a spin-3/2 hole, 71 and 72 
are the Luttinger parameters. Note that here the rela- 
tion between Vz and B differs from Eq. ([T]) by a factor 
1/2 because we use the spin-3/2 matrix J^, instead of the 
Pauli matrix which was used for Eq. ([1]) [13] (thus the 
required B for the same Vz is smaller by a factor 1/2). 

To simplify the calculations we assume a rectangular 
cross-section of the nanowire with the widths Ly and 
Lz. The strong confinement along the y and z directions 
makes the energy levels quantized on these axes. To il- 
lustrate the emergence of the Majorana fermions, we first 
consider a single band, i.e., the lowest energy state along 
the y and z directions. Using the ground state wave- 
function ^(7/, z) along the y and z directions, the original 
Hamiltonian (|2]) in 3D can be projected to an effective 
ID form 



H^{x) = (7i/2 + 572/4- 72Jj)9^+^S^'/^' 

Wj2Jz/Ll + iaJyd^ + VzJz - ^ - (3) 



where ^ = 5727r^ [Ly ^ + ^) /4, and fi = fi ^ 
7i7r^ {Ly'^ -\- 1/^^) /2 is the shifted chemical potential due 
to the confinement. 

In Fig. [Ud, we plot the energy spectrum ss of the 
Hamiltonian (|3]) for holes in an ID geometry. There are 
two heavy hole and two light hole bands. If fi lies in 



the shaded region, it intercepts only one Fermi surface. 
An odd number of Fermi surfaces implies a breakdown 
of the fermion doubling theorem (due, in this case, to 
the Zeeman splitting), which yields, in the presence of a 
superconducting pair potential, the required TS state for 
Majorana fermions [10|. The position of the shaded re- 
gion (the TS state) is determined by the Hamiltonian ([3]) 
at A: = that depends on the parameters 71, 72, Ly^ Lz, 
while the width of the shaded region is determined by Vz . 
The large effective mass and strong spin-orbit coupling 
of the holes lead to a high density n = f^^ dk/2TT ~ 10^ 
cm~^ of holes in the TS state. Such a hi gh density is rou- 
tinely realized in nanowire experiments |20M22] . The high 
carrier density provides a strong screening of the disor- 
der potential, suppressing the spatial chemical potential 
fluctuations and disorder effects [23]. 

The superconducting pair potential can be induced 
in the hole-doped nanowire through the proximity con- 
tact with an 5- wave superconductor (Fig.[TJi), as demon- 
strated in experiments [22^. This yields the Hamiltonian, 



Hsc = Y. I d^'r^smj (r) V'L.^L, + H.c, 

^ — ^mj J J J 



(4) 



where V^J^^ are the creation operators for holes with the 
angular momentum mj = ^, f and A^^^ (r) is the prox- 
imity induced pair potential. The form of the pairing 
Hamiltonian is dictated by the fact that A^^^ couples 
particles with mj with particles with —mj, and should 
be determined through the microscopic theory of the 
proximity effect [26] . 

Taking account of the spin-3/2 and the particle- 
hole degrees of freedom in the superconductor, the 
Bogoliubov-de-Gennes (BdG) Hamiltonian can be writ- 
ten as an 8 X 8 matrix 



^BdG = 



A*(x) 



-Tti^*(x)T 



(5) 



in the Nambu spinor basis l>(x) = ^?/;(x), T?/;'''(x)^ 

with T = i{l2 ^ cFx)Ty {I2 is the 2x2 unit 
matrix), A5(x) = diag(A53/2, A51/2, A51/2, A53/2) , and 
^p{x) = {tjj3{x),2Pi{x),2p_i{x),ip_3{x)). 

Parameter space for the topological state: In a 
ID nanowire, the parameter regime for the Majorana 
fermions (in ID the Majorana fermions are localized at 
the two end points) can be determined by the topological 
index X [23, [H defined as. 



M = sgn [Pf {r(0)}] sgn [Pf {r(7r/a)}] . 



(6) 



Here Pf represents the Pfaffian of the anti-symmetric 
matrix T{k) = —iHBdG{k){^y T), c^y is the Pauli ma- 
trix, HBdcik) is the corresponding BdG Hamiltonian in 
the momentum space {—idx k), and a is the lattice 
constant. M. = —1 (+1) corresponds to the topologi- 
cally nontrivial (trivial) states with (without) Majorana 
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FIG. 2: (a) Plot of the phase boundary = 0) between 
topological and non-topological states for different /x, Vz and 
Tj. Asi/2 = 1 meV, ^2 = 55.4 meV (corresponding to Ly = 10 

nm), /3i = (1 - ry)/(2v1^)^/3|-ryA2^/2- The color region 
represents the TS state with M— —1. Different colors cor- 
respond to different critical Zeeman fields, (b) Plot of the 
minigap Eg with respect to A. /i = — 50.9 meV, Vz — 1.5 
meV, Ly ^ 9 nm and 15 nm. 



fermions. Using the fact that T{k) is an anti-symmetric 
matrix that can be diagonaHzed by a lower triangular 
matrix [29^, we find 

Pf {r(fc)} = Pf {AsT} Pf {TA5 + Hf (k) TAg'H, (k)} 

(7) 

through a straightforward calculation. 

In the continuous limit /c = 7r/a ^ 00, the terms 
in the single particle Hamiltonian Hi (k) dominates and 
all other terms in T{k) can be neglected. In this case, it 
can be shown that sgn[Pf{r(A:)}] = sgn[det {Hi (k))] = 1, 
therefore M. is solely determined by the sign of Pf{r(0)}. 
Pf{r(0)} can be derived analytically from Eq. (O, yield- 
ing A4 = sgn[J^], where 



^=/o-/iV;' + 9^^/16, 



(8) 



AS3/2A51/2 

- /3l(A,3/2 

+ 9A? 



^s3/2 



6/?| 



/i 

/4, 



\/37r272i-V2. Since 



/o = + 

(As3/2 - A51/2)// ^ 
10^2 + ml + l&Plfi -r ^-,1/2 
/?1 = 7r272(i:j2_£-2/2), ^2 

M changes sign when changes sign, the phase bound- 
ary between the topologically trivial and nontrivial 
states can be determined by setting T = {). 

When the ratio r] = A53/2/A51/2 > Al^^^/f3^, there 
exists a magic value of the chemical potential fto = 

/o = 0, where /3i = {1 - rj) fto/ {I ^ v)- For /o = 0, 
< and Major ana fermions exist even for a vanish- 
ingly small Vz. Clearly, r] > A^//^! requires that Ag3/2 
and Asi/2 have the same sign and /^l > A1A2. Be- 
cause of the strong confinement, the second condition 
can be easily satisfied. Therefore the threshold Zeeman 
field V^ for the TS state vanishes when As3/2 and Asi/2 
are of the same sign, independent of their relative mag- 
nitudes. When As3/2 and Asi/2 have different signs, V^ 



77/31 — ^33/2 heavy hole band such that 



becomes nonzero, but is still much smaller than that for 
the electron-doped semiconductors. Therefore the rel- 
ative signs of the pair potentials should not matter in 
realistic experiments because a reasonable Vz is always 
needed to generate a sizable chemical potential region 
for the TS state. In Fig. 2, we plot the boundary = 
between topologically trivial and non-trivial states for 
different /i, V^ and r]. The TS states for a fixed V^ are 
embraced by two lines with the same color for the cor- 
responding V^. For instance, for 77 = 0.15 and V^ = 1 
meV, Fig. 2 shows that the TS state exists in the region 
111 < ji < II2 with jii = —103.3 meV and 112 = —100.8 
meV. Clearly, the Vz = line (the center blue line) ex- 
ists for 7^ > 0, but vanishes for 7^ < 0. When 7^ ^ 0+, 
j2o —00. 

The vanishingly small V^ for the TS state at 7^ > 
may also be understood by projecting the four band 
Luttinger model in Eq. (|3]) to an effective two heavy 
hole band model because of the large energy splitting 
between the heavy and light hole bands. The resulting 
two band model is generally very complicated for finite 
kx and Vz. However, because the Pfaffian is determined 
by the Hamiltonian at /c = and we are interested in 
the TS state with vanishingly small Vz^ we can do the 
band projection at 14 = and around kx = 0, leading 
to an effective pairing A^^^ = (A53/2 — tzAgi/2) /i^ and 
an effective chemical potential peff = P — + /^l , 

with f<i = (^\/l3l/l3l + 1 - pi/pQ^ . Here we have ne- 
glected non-diagonal term of the pairing because /^l ^ 
A53/2A51/2. When 7? > 0, A^^^ may vanish by choosing 
/3i//32 = (1 — ^) /2y^, therefore the critical Zeeman field 
V^ also vanishes when /ig// = 0. While when 7^ < 0, 
A^^^ is always finite and there is a minimum V^ based 
on Eq. ([1]) for the two-band model. Note that the zero 
A^^^ at /caj = and Vz = does not imply the zero A at 
a finite Vz and kx. For a large kx^ the coupling between 
heavy and light holes becomes important and the mini- 
gap is finite for a large Vz even at the magic /io (see Fig. 
3b). 

Henceforth, we consider two representative cases to 
further illustrate our results: (I) A53/2 = Asi/2 = A 
and (II) A53/2 = — Asi/2 = A. For the a reasonable 
large Vz^ both pairings yield similar observable signals. 
The corresponding topological region is plotted in Fig. 
3a. In the case (I), the TS states exist in the region 
\P~ Po\ ^ 1^^1/2 for both positive and negative Vz. The 
transition at 14 = (at which the superconductor is 
gapless and non-topological) is a quantum transition at 
which neither the symmetry of the system nor the topo- 
logical properties change as a function of Vz. Note that 
the system crosses the phase transition boundary lines 
twice when Vz changes from —00 to -|-oo for a fixed /i, 
except at the magic value po- At the phase boundary, 
the quasiparticle energy gap closes and the superconduc- 
tor becomes gapless. At the magic /io, the two lines of 
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FIG. 3: The parameter regime for the existence of Majorana 
fermions. (a,b) and (c,d) correspond to the single and two 
band models. Ly ^ 10 nm and Lz ~ 14 nm for (b) and case 
(I) in (a). Ly ^ 9 nm and Lz ~ 15 nm for (d) and case (I) 
in (c). Ly ^ 14.2 nm, Lz ~ 9.5 nm for case (II) in (a) and 
(c). The widths of the nanowire are chosen to obtain a small 
critical Zeeman field for each case. (a,c) are obtained from 
the topological index M = —1 (the filled regions). (c,d) are 
obtained from solving the BdG equation with Vz = 1.5 meV. 
A = 1 meV. The dashed lines in (b,d) give the minigap. The 
parameters 71, 72, a are the same as that in Fig. lb. 



the phase boundary merge at = 0, therefore the su- 
perconductor is gapped and topological for all Vz except 
at Vz = 0. 

In the case (II), the threshold V^ ~ pA with p = 
2 [1 + 2 (1 + (see Figs. 3a and 3c). For 

instance, for hole-doped In As nanowires (with a typical 
gl^ = 35 [30]) with Nb as the adjacent superconductor 
(A 1 meV), the required magnetic field B is ^ 0.35 T, 
which is about 1/3 of the corresponding B for electron- 
doped InAs nanowires (with the same material Nb as the 
adjacent superconductor) [17], and can be easily realized 
with a bar magnet without affecting the superconduct- 
ing pairing. Therefore, we find that irrespective of the 
relative sign of and Asi/2 (which can only be de- 

termined from a more microscopic calculation [26]) the 
threshold B for the TS state in the hole-doped case is 
much smaller than that in the electron-doped case. The 
fact that the Majorana fermions can be observed even 
with a small magnetic field opens the possibility of using 
a wide range of semiconductor materials with only small 
gl^ factors [21, 22] and is one big advantage of using hole- 
doped semiconductors. 

To further confirm the existence of the Majorana 
fermions in the above parameter regime, we also numer- 
ically solve the BdG equations (|5|) and obtain the energy 



spectrum and eigenstates. The Majorana fermion corre- 
sponds to a zero energy eigenvalue in the BdG spectrum. 
Henceforth, we present our results only for the case (I), 
but have confirmed that the results for the case (II) are 
similar. In Fig. [3]3, we plot the ground and the first 
excited state energies. The ground state energy becomes 
zero in the region /ic2 < /i < /ici- In the same param- 
eter region, the topological index A4 = —1 (Fig. [3ti), 
agreeing with the numerical results. The solution of the 
BdG equation also yields the minimum energy gap (mini- 
gap) above the zero energy states. At this gap and above 
there are other, finite-energy, states localized at the end 
points of the wire. The minigap therefore provides the 
protection for the Majorana states from finite tempera- 
ture thermal effects. We see that the minigap is of the 
same order of the pairing gap A 1 meV, which means 
that the Majorana fermion physics can be accessed at the 
experimentally accessible temperatures T < 10 K. 

Effects of multiple confinement bands: In a realistic 
experiment, multiple confinement energy bands along 
the y_ and z directions need be taken into account 



[18|, |3l| because they are mixed with each other by 
the large spin-orbit coupling. Considering the lowest n 
confinement bands, the BdG Hamiltonian can be written 
as an 8n x 8n matrix similar as Eq. ([5]) with the matrix 
elements replaced with the corresponding multiband 
forms. Specifically, Hi (x) is replaced with Hn{x) = 

where (j)i{ii^z) is the wavefunction on the i- 
th band in the yz plane. T is replaced with 
= In ^ T. The Nambu spinor basis becomes 
^(x) = {M^)r-- ^Mx),Tyl{x),--- ,TV^i(x))^ with 
2pi{x) = {ip3.{x),2pi^{x),ip_i^{x),ip_3.{x))^ as the 
hole annihilation operator on the i-th band. We 
also assume there is no superconducting pairing be- 
tween holes at different confinement bands, therefore 
Asn = In ^ As. Following a similar procedure as that 
for the single band, we derive the topological index Ai = 
sgn[Pf{As„T„}Pf{T„As„ + (fc) T^A^^H^ (k)}] 
for the multiband model. 

Here we consider only the lowest two relevant energy 
bands {ly = 1 and 2, = 1). In Fig. 3c, we plot the 
parameter regime for the TS state when the multiple 
confinement induced bands are included. The parame- 
ter regime for the TS state now has some quantitative 
difference from that in the one-confinement-band case. 
However the basic conclusion remains the same, that is, 
the Majorana fermions can be realized even with very 
small Zeeman fields. In Fig. 3d, we plot the energies 
of the ground and first excited states by solving the rele- 
vant BdG equations in the multiband model. The figures 
confirm that the two methods - topological index and 
numerical solutions of BdG equations - yield the same 
results. In the multiband model, the minigap is slightly 
reduced, but still at the same order of A. 

In practice, A may depend on the material as well 
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as the interface between the superconductor and the 
nanowire. Although the parameter regime for the emer- 
gence of Major ana fermions does not change much as a 
function of A, the minigap has a strong dependence on 
A. In Fig. 2b, we plot the minigap with respect to A 
in the multiband model and find that Eg ^ A (simi- 
lar as the electron-doped case p^), instead of A^ as in a 
regular 5- wave or a chiral-p wave superconductor. There- 
fore the minigap is rather large, which ensures thermal 
robustness of the Majorana fermions. 

Conclusion: The list of systems capable of support- 
ing a non-Abelian TS state now includes the filling 
factor u = 5/2 FQH state, chiral-p wave supercon- 
ductor s/superfluids, topological insulators, and electron- 
doped semiconductors. To this list we have added a 
new system, a hole-doped nanowire, as a possible non- 
Abelian platform. Although the roster is growing, ours 
is not an ordinary addition. As we have shown here in 
detail, the requirements (carrier density, magnetic field, 
^-factor, etc.) for the TS state in a hole-doped nanowire 
are already accessible in experiments. Thus this system 
can be a potential breakthrough facilitating solid-state 
demonstration of Majorana fermions as well as realiza- 
tion of TQC using a nanowire network. 
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